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Abstract: A dual component made of non-relativistic particles and a scalar field, ex- 
changing energy, naturally falls onto an attractor solution, making them a (sub)dominant 
part of the cosmic energy during the radiation dominated era, provided that the constant 
/3, measuring the coupling, is strong enough. The density parameters of both components 
are then constant, as they expand as a -4 . If the field energy is then prevalently kinetic, as 
is expected, its energy is exactly half of the pressureless component; the dual component 
as a whole, then, has a density parameter VL c d = 3/4/3 2 (e.g., for j3 ~ 2.5, fl c d ~ 0.1, in ac- o 
cordance with Dark Radiation expectations). The stationary evolution can only be broken 
by the rising of other component (s), expanding as a~ 3 . In a realistic scenario, this happens 
when z ~ 3-5 x 10 3 . When such extra component(s) become(s) dominant, the densities 
of the dual components also rise above radiation. The scalar field behavior can be easily 
tuned to fit Dark Energy data, while the coupled DM density parameter becomes 0(1CP 3 ). 
This model however requires that, at present, two different DM components exist. The 
one responsible for the break of the stationary regime could be made, e.g., by thermally 
distributed particles with mass even 3> 1-2 keV (or non-thermal particles with analogous 
average speed) so accounting for the size of observed galactic cores; in fact, a fair amount 
of small scale objects is however produced by fluctuation re-generated by the coupled DM 
component, in spite of its small density parameter, after the warm component has become 
non-relativistic. 
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numerical. 



© SISSA/ISAS 2013 



http://jhep.sissa.it/J0URNAL/JHEP3.tar.gz 



Contents 

[l]. Introduction |l| 

g. CDM-DE coupling | 

|3|. Coupled DE in the radiative era |5] 

g. Stability § 

[5]. A Lagrangian approach ||| 

|6]. Exit from the radiative regime 10 

[7]. Astrophysical context: flat galaxy cores and small halo deficit 13 

pL Discussion 16 



1. Introduction 

Cosmological models supported by data are affected by a number of fine tuning and co- 
incidence paradoxes. Among them, we have the problems concerning Dark Energy (DE), 
namely if its state parameter w = —1 : (i) Why are we living in the only era when DE 
and matter have comparable densities? (ii) Why density inhomogeneities are normalized 
so that DE allows fluctuations to reach the non-linear regime, then stopping any further 
density evolution? 

Models trying to ease these paradoxes were proposed in the last few years, often based 
on the idea that DE is a scalar field 0, @, §, possibly coupled to Dark Matter (DM) |§. 
Current data, however, do not provide sufficient discrimination, so that no such model 
really appears statistically favored in respect to ACDM (see, e.g., ||). Accordingly, rather 
than suggesting models, recent work has focused on planning measures allowing us to dis- 
criminate between a DE state equation w{a) = — 1 and other behaviors, as the forthcoming 
Euclid mission 1 |j] 

In this paper, in a sense, we partially go back to the older approach, by suggesting 
a class of cosmological models, based on the assumption that DE is a scalar field <j). In 
describing them, however, no peculiar self-interaction potential V((p) is selected, as the 
features we point out are (almost) independent from it, so that changing V((p) has a 
modest impact on our findings. Accordingly, while we expect that the shape of V((j)) can 
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be hardly discriminated from ACDM through CMB and related data fits (see, e.g., @), 
we suggest specific predictions possibly discriminating this class of models from different 
cosmologies. 

The basic point made in this paper is that a large fraction (typically 1/3) of Dark 
Radiation (DR) is a scalar field <fi, later turning into Dark Energy. Big-Bang Nucleosyn- 
thesis (BBN) is consistent with observable nuclide abundances in the presence of an extra 
radiative component consistent with AN jf < 1.26 neutrino species ||, and this constraint 
can be furtherly softened in the presence of primeval lepton anti-lepton asymmetry. CMB 
data are slightly more constraining, as the best fit of WMAP and related data yields 
AN ff = 0.85 ± 0.62 (2 standard deviations) [[R]]. A large deal of work has therefore 
been dedicated to the nature of DR, also outlining that the amount of DR could be quite 
different in BBN and CMB data Q. 

The price to pay for DR turning into DE is that two kinds of Dark Matter (DM) must 
exist, that we shall dub DM cou and DM unc . The former one (DM cou ) is coupled to the 
scalar field 4>. In the radiative era it is the rest of DR; in the Newtonian limit, it has 
ordinary gravitational interactions with any other cosmic components, but its fluctuations 
feel a much stronger self gravity, while its dynamics has further specific modifications |l2]] . 
Its density parameter, in the present epoch, shall lay in the per mil range, but its peculiar 
gravitational behavior can give it an important role in shaping a large deal of today's 
observables. The latter DM kind (DM unc ), expected to have an ordinary gravitational 
behavior, is then needed to break a primeval stationary condition. In principle, quite a few 
discriminatory predictions may follow from the presence of the DM cou component, coupled 
to DE. 

Although the mechanism turning DR into DE is the basic point of this work, we shall 
devote a specific Section to discuss possible scenarios allowed by the presence of 2 DM 
kinds. Rather than a complication, this appears as an "opportunity". In particular, we 
shall discuss the case of DM unc being warm; the simultaneous presence of DM cou might 
then ease a number of cosmological problems. Among them, the small halo deficit and the 
size of the plateau in galaxy cores, when comparing observations to simulations. 

All that is obtained by pushing to extreme consequences the idea of DE-DM coupling, 
as previously suggested in The coupling intensity considered is however much greater 
than any previous analysis, and consistency with data can be recovered thanks to the fact 
that DM cou is a minor component of DM. 

The plan of the paper is as follows: In the next Section we shall discuss how to treat 
coupled DM and DE if the DE state equation w(a) is assigned, while the DE self-interaction 
potential V((j)) is unknown. The equation of motion to be solved is then simpler than the 
usual Klein-Gordon equation, as the problem becomes first order, the real unknown being 
4>\ = 4>, while we do not need to know (p. In Section 3, we then consider the 4>\ equation and 
the equation ruling DM cou density in the radiative era, finding a self-consistent solution, 
which allows constant density parameters for DM cou and 0-field, provided that the coupling 
strength is large enough. In Section 4 we verify this solution to be an attractor, and 
that we converge on it if starting from any initial condition. In Section 5 we briefly 
discuss the form of the effective Lagrangian yielding the equation of motion. Section 6 is 



then devoted to debating what happens if another non-relativistic component finalizes the 
radiative expansion. In Section 7 we then show that, rather than a complication, the two 
DM components are an opportunity, in particular if we assume DM unc to be WDM. A 
Discussion Section concludes the paper. 



2. CDM-DE coupling 

The possibility that CDM and DE are coupled have been considered by several authors 

[l3|, [14]]. As a matter of fact, while the stress-energy tensors of CDM and DE, 
and ^ respectively, surely fulfill the pseudo-conservation equation 



T (cC + ^w = o, (2.1) 

there is no direct evidence that the two equations 

T ( C w = °' r <«C = ( 2 - 2 ) 

are separately satisfied. The r.h.s. of the above equations can then be replaced, in a 
covariant way, by a term yielding a leakage of energy from DE to CDM or viceversa. 
Here we shall consider the option yielding 

4> + 2-<p + a 2 VM) = Cp c a 2 , p c + 3-p c = -C# c . (2.3) 
a r a 

Here p c is the density of CDM, eft is a scalar field self-interacting through the potential 
V{4>) and accounting for DE, while 
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accounts for energy transfer from CDM to DE. Taking a constant B is an extra assumption 
we shall make here for simplicity. Here we also assume a FRW metric reading 

ds 2 = a 2 (T)(dT 2 -df) , (2.5) 

so that a is the scale factor, dl is the spatial line element, while differentiations are made 
in respect to the conformal time r. 

The option (2.4) has a peculiar significance, as it allows models where DE is always 
a significant component of the Universe, being fed energy by the CDM component which, 
accordingly, dilutes more rapidly than a~ 3 . 

Let us then reconsider the eqs. ( |2.3| ) when the potential V{<f>) is unknown, while we 
know that the DE state equation is a suitable w(a). It must however be 

<P 2 /2a 2 - V . <j> 2 1-w 
w = , i.e. V = 7T- , (2.6) 

and a 2 V^, in eq. (|2.3| ), can be directly evaluated by considering how the two factors in the 
last expression depend on r, then associating such dependence with the r dependence of </>. 
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d_l_ 

dtft 2 

d ( 1 1-w 



da \a 2 1 + w 



d 6 d6 dr 

' y — x 

ddi 2 dr dd> 



1 — w dw 
+ 



(2.7) 



so that 



Let then 



a 2 Vi 



1 — w a 



1 + w a 
1 



l + w da (1 + wf 



1 — w dw 
+ 



l + w da (1 + w) 2 
a dw 



W = - I 1 + 3w 

2 V l + w da 



in order that the system of equations (^^) becomes 



h + W^i = ^Cp c a" 
a 2 



p c + 3-p c = -Ccfripc . 
a 



(2.9) 
(2.10) 

(2.11) 



Here we set (f>\ = <f> to outline that the former equation has become first order. 

As a matter of fact, it seems more realistic that data allows us to know w(a), rather 
than the potential V(<p). Should we know w(a) and wish to interpret the w(a) dependence 



as due to the evolution of a scalar field, we need to integrate just eqs. (|2.1l| ), together with 
the Friedman equation. The scalar field contribution to its source term is then the total 
energy density for DE, 

(2.12) 



Pd 



i 2 (l + w) 



Appar- 



as also the potential contribution is derived from <j)\ and w(a), thought eq. (^ 
ently, therefore, we need not recovering the un-differentiated <f> behavior. 

This however assumes that we know DE and DM cou to be coupled, that the coupling 
is constant, and the coupling constant has a specific value C = 4(7r/3) 1//2 /3/m p . 

It is premature to discuss here observational strategies. A natural start point, however, 
amounts to assuming no coupling. The apparent DE state parameter, measured from from 
the expansion rate, would then be 



with 



w eff (a) = w(a)/[l + £(a)] 



£(a) = \g(</>)/g(<h) - 1] * /W(/W) 



(2.13) 



(2.14) 



Here g{4>) oc p c a 3 tells us the deviation of p c from the ordinary a -3 scaling; the suffix o 
refers anywhere to today's quantities (see || for a detailed discussion). 

Notice that, in order to pass from w(a) to w e jf(a) or viceversa, we then need to know 
4>, besides of <j>\. However, to recover (f>, we do not need integrating an equation, but just 
a known function. 

We shall not delve here into a possible more refined analysis, seeking the family of 
wp(a) behaviors as the assumed coupling j3 varies. 



o 
o 

N) 
o 
i — 1 

rO 

O 
O 

o 



- 4 - 



Notice also that eq. ( 2.12 ) shows that p<i is positive definite only if w > — 1: to explore 
the w < — 1 domain one needs the (f> field to have a suitable anomalous kinetic energy 
expression. 

These equations, therefore, require no specific potential shape to be assigned, no back- 
ground expansion regime to be assumed, while the very w(a) behavior is generic, although 
it must be w(a) > —1. 

3. Coupled DE in the radiative era 

Let us now consider the system ( |2.11| ) when the background expansion is supposed to be 
radiative. We shall make the further assumption that 

W = ±(l + 3w), (3.1) 

with w = const., as is reasonable when a — > 0. This assumption is however unessential and 
only allows us to simplify the analytical treatment. 

Let us also remind that, in the radiative era, it is a oc r, so that a/ a = 1/r. Accordingly, 
from Friedmann equations we obtain that 

*~ paV = l, (3.2) 



3m^ 



p being the background energy density and m p the Planck mass. The latter eq. ( |2.11| ) has 
then the formal integral 

Pc = Pi,c(—) exp ( -C [ drfa) , (3.3) 



Ti being a reference time when CDM density is p c j and the scale factor is a% = a(rj). 
If this expression for p c is then replaced in the former eq. ( 2.11| ), we have a first order 



transcendental differential equation whose unknown is <j)\. It seems hard to find a generic 
analytic integral of this equation. 

There is however a peculiar case, allowing integration. Let us make the ansatz that 

01 = . (3.4) 
r 

Taking eq. (2.4) into account, we have then that 

-C [ dr4> x = ln(-) ab , (3.5) 



so that 



T 



3+ab 



Pc = Pi, c y—j (3.6) 
and, by replacing the expression ( |3.4f ) in the former eq. ( [2.11 ), we obtain 



m„ 1 + w b fa T \ 3 + ab 

- = — ^ (") • < 37 » 



In order that the two sides scale with a in the same way, it must then be ab = 1 and the 
DM density shall scale with a -4 . The fact that p c dilutes more rapidly than oc a~ 3 does 
not come as a surprise, as there is a continuous leakage of energy from it to the (j) field. 
The fact that it dilutes exactly as a -4 , instead, is a consequence of the ansatz Q3.4| ). 
Equation ( |3.7|) can then be put in the form 

1 W - 1 8tt 2 2 

■p c ar =S1 C , (3.8) 



/3 2 1 + w 3m| 

owing to eq. (fT^). Here S7 C = p c /p is the (constant) density parameter of DM, during 
radiation era. In order that the ansatz ( |3.4| ) is allowed, O c ought to have the value given 
by this equation. 

Also the energy density pd of the DE field (j) scales with a -4 . In fact, owing to eq. ( 2.12| ), 



a 2 ml 1 

Pd = -227-; — > ( 3 - 9 ) 

and, using again eq. fl3,2| ), we obtain the constant density parameter of DE 
showing also that 

^£ = 2(W - 1) = 3w- 1 . (3.11) 

Accordingly, the whole framework is consistent only if w > 1/3. This means that the energy 
density of the kinetic part of the <j) field should be dominant in respect to the potential 
part. In the specific case w — 1, holding for (f>\ 2a 2 V, we therefore expect that it is 
constantly 

n c ~ 2 U d (3.12) 

during such expansion regime. 

Altogether, these computations show that, during a radiative expansion, we can have 
two coupled components also expanding oc a~ 4 although their state equations are w c = 
and w > 1/3. The former component is made of non-relativistic particles, the latter is a 
self-interacting scalar field. A possible option is that the individual CDM particle masses 
decrease in time, as a consequence of their interaction with <j>. The two dark components 
interact with a strength measured by the dimensionless parameter /3. 

Owing to their behavior, these two components do not modify the radiative charac- 
ter of the expansion. However, we might prefer that ordinary radiation is the dominant 
component during this period. This requires a strong coupling between the components 
/3 S> 1, as it should however be w < 1 . In turn, the ratio between the density of DM and 
DE is 0(1). 

In the specific case w ~ 1, however reasonable in the very early Universe, we then have 
n d (] 2 = 1/4 and 

(n c + n d )p 2 = 3/4 . (3.13) 



a 
+ 

o 







-0.5 - 



-1 



-1.5 - 



1 


1 1 1 
■a 


U 1 1 1 1 | 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 


- 


allowi 




w= + 1 

VV 1 ± _ 


- 


not 




- 




—> 
s 








0) 

c 








o 
a 








£ 
o 








o 








"3 








3 
T3 






1 


1 1 1 


1 1 1 1 1 1 1 1 


i i 1 i i i i 1 i i i i 



1 2 3 4 5 

Figure 1: Density parameter of the dual radiative component vs. the coupling parameter (3 



In an early epoch it is fair to assume a total density parameter 0^ = 1. Then, requiring 
Q c + Qd < 1 yields 

P > V3/2 = 0.866 . (3.14) 

A solution with w — 1 and /3 2 ~ 3/4 , although self-consistent, seems unreasonable. In 
fact, then £l c + ~ 1 and the ordinary radiation component should vanish. 

Solution with (5 2 < 3/4 require u; > 1 to allow constant 0, c ,d- If it? > 1 is excluded, when 
(3 2 < 3/4 there exist no solution with constant tt c ,d, i-e., the CDM and 4> field contributions 
to the overall density become increasingly small when a tends to zero. 

In Figure |l| we plot the density parameter of the dual (DM+DE) component vs. as- 
sumed j3 values. 

The dual component gives place to a natural form of DR. In general, we can gauge its 
significance through the number of extra neutrino species 
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(here pdr is the dual-component density and T is photon temperature). It is then 

(8/7)(ll/4) 4 / 3 +3 7.4032 



(3.15) 



AN, 



off 



(4/3)02 



1 



(4/3)/?2 



1 



(3.16) 



and P = 1.5 (2.5) yields ~ 3.7 (1.1) extra species. 

Let us recall that BBN prescribes AN Q ff < 1.26, in standard theories, or AN Q ff < 2.56 
when a non-vanishing chemical potential is allowed for neutrinos ||. The overall density 
of radiative components also sets the equality redshift. CMB and related data analysis 
then allow us to fix the equality, so requiring AN ff ~ 0.85 ± 0.62 or similar figures [pTOf] , 
although the choice of priors risks to affect the final estimate fi~5||. All above limits are at 
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Figure 2: Stability of solutions. We start from modified initial values, showing that we soon 
reconverge on the "attractor" solution. More specifically, we shifted by ±50%, respectively: p c , 
<fii, both of them; the last option is equivalent to setting a "modified" j3 value. The solid lines are 
the "attractor" solution. Dashed and dotted lines show the gradual recovery of it. In these Figures 
13 = 2.5 . 

95 % confidence level. Let us also recall that models can be easily built, where the DR 
density at equality exceeds the one at BBN. The smallest coupling consistent with all above 
limits is however around (3 = 2. Through this paper, only the case f3 = 2.5 will be however 
considered. 



4. Stability 



Solutions with w = 1 and 

fl c = 2Vt d = 1/2/3 2 



(4.1) 



are however stable. If the initial values of S7 C , 0, or f3 do not fulfil the above relation, and 



the expansion regime is radiative, the condition (4.1) is soon restored. 

In Figure || we show this in 3 cases: (i) If we set an initial value of p c in excess by 50 %. 
(ii) If we set an initial value of 0i in excess by 50 %. (iii) If both shifts are simultaneously 
performed; this is equivalent to having a "wrong" initial (5. 

The results in Figure Q are obtained by numerically integrating the set of differential 
equations fl2.3| ) plus the Friedmann equation 

(a/a) 2 = 8vrp/3m 2 , (4.2) 

p being the total background density. All quantities are expressed in MeV; in particular 
the units for the conformal time r (not shown) are MeV -1 . 

The above numerical output can be analytically understood. Let us consider, e.g., that 
the density of the DE component has a value different from what eq. (3.10) requires, i.e. 
that 

01 = 0l,o + 8 an d 0i,o = arrip/T (4-3) 

with 5 positive or negative. Assuming that p c is unmodified, eq. ( 2.11 ) tell us that 5 must 
fulfill the equation 

8 + (1 + e)8/r = (4.4) 
with e = W — 1 > 1 (for w > 1/3). Then, 5 oc r~( 1+e ) and the ratio 

|<5|/0i oc r~ e (4.5) 

necessarily decreases with time. As a matter of fact, however, 0i appears also in the 
differential equation setting p c , while the impact on Friedmann equation is small if we 
assume a large (3 yielding small Q Ct d (plots are for f3 = 2.5). A greater 0i then yields a p c 
decrease faster than a -4 . In turn this corresponds to a decreased energy leakage towards 



the scalar field, so that \5\ decline is accelerated, in respect to (f4.5|) . Figure |2j however 
shows that some bounces occur before the "attractor" solution is recovered. 

Analogous arguments can be put forward for the cases (ii) and (iii). For the sake 
of completeness, let us also outline that solutions are tendentially stable also for constant 
w < 1. However, if initial conditions violate eqs. (3.10)-(3.11), the recovery of the attractor 
solution takes an increasingly longer time as w is farther from unity. 

5. A Lagrangian approach 



The function g(0) needed to define £ in eq. ( 2.13 ) also enters in the Lagrangian coupling 



between the scalar field and a supposed spinor field ip, yielding DM cou . 

According to ||, the general expression of the effective interaction Lagrangian reads 

C = , (5.1) 

p, being a factor with the dimensions of a mass. Owing to eq. (|3.3D, then, it should be 



g(0) = exp -b / dT<pi/m p (5.2) 




Figure 3: Density evolutions when a non-relativistic component becomes dominant at z ~ 5 x 10 3 . 
Both DM C£m and DE densities then overcome radiation, but keep steadily below DM mlc . We 
assumed here (3 = 2.5, so that J7 C + ild = 0.12 in the radiative era; this corresponds to AN a fj = 
1.009, as Dark Radiation, i.e. to ~1 extra neutrino species. 

so that, owing to the ansatz ( |3.4D , 

C = u—tpip = b-^^tp = T^yipip , (5.3) 
r 

with a constant T whose dimensions are to -1 , or 

C = fL-fyjj . (5.4) 
m p 

This makes clear that the very interaction Lagrangian also displays the role of variable- 
mass term, with a mass oc r _1 or oc T. 

6. Exit from the radiative regime 

The picture changes if the expansion is no longer radiative. In the physical world, there 
will be baryons, at least, whose density overcomes the density of the radiative components 
slightly below z ~ 10 3 . If their abundance is consistent with BBN, they are surely not 
enough to produce an overall picture possibly close to observations. We shall therefore 
assume that another non-relativistic DM component exists so that its density summed to 
baryons matches the density of the radiative component at z = 5 x 10 3 . We have therefore 



two DM components: DM cou coupled to the field </>; DM unc uncoupled; for the sake of 
simplicity, here baryons are included in DM mc . 

In Figure ^| we show how DM cou and the DE field no longer scale as radiation, when the 
overall expansion ceases to be radiative. We assume that (3 = 2.5, so that AN ff ~ 1 , and 
that the equation of state of DE keeps w = +1 until today. Then, (j) does not contribute 
much to today's overall density while, by decreasing /3, we could have DM unc approaching 
the present baryon density, at most. Below, we shall further comment on DM nrac observable 
effects. 

The situation changes radically if we include the option that the (j) field abandons the 
kinetic regime. Of course, this is necessary if we wish to identify it with DE, whose today's 
state equations approaches w ~ — 1. For most potentials considered in the literature, when 
4> overcomes a suitable level, the potential energy becomes indeed dominant. The moment 
when the transition occurs depends on the potential assumed as well as on the initial 
contribution of the two coupled components to the cosmic budget. 

As an example, we report here the expected behavior of w(z) when the potential is 
SUGRA Q or RP |E|. More details on these plots are given in Q. The Figures f| show 
the shift of w from +1 to ~ —1, with some potentials studied in the literature. Values 
p < V3/2 are however taken in these Figures, while DE is coupled to all DM. 

Our aim here, however, just amounts to showing that DE can achieve a density con- 
sistent with observations while DM con , on the contrary, yields a negligible contribution to 
the cosmic budget. This unavoidably requires that the energy density of the <fi field turns 
from (prevalently) kinetic to (prevalently) potential. 

The only caution to be taken is avoiding a too fast w decrease, as the expression of W 
has a contribution from dw/da, which may become dangerously high. 

In this work we take the following class of interpolatory functions 



1 + w = (1 + uu + ) exp 



for a < a\ 



l + w = l±^±(^l) fora>ai (6.1) 



a 



yielding 

a dw ( a 



1 + w da \a\ 
a dw 



< e for a < a\ 



-e for a > ai (6.2) 

1 + w da 

so that the W correction, in respect to 1 + 3w is e, at most. Here w+ is the DE equation 
of state at large z. In principle, e is to be fixed so that the DE state equation is a suitable 
W- at z = 0. Here we used e = 1.9. In Figure [| we show the resulting w(a) behavior, when 
l + zx = l/ai = 12 . 

Notice that the w(a) behavior shown in the Figure is not so far from those obtained 
from some assigned potentials. In particular, the large value of /3 is balanced by the low 
DM cou density which, however, is not arbitrary, as the initial T)M cou density parameter is 
dictated by j3 itself. 
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Figure 4: Scale dependence of the DE state parameter in SUGRA and RP models, in the presence 
and absence of energy transfer from a large CDM component. The main parameters of the models 
are n b = 0.046, fl c = 0.209, H = 73 km/s/Mpc, while A = 0.1 GeV for SUGRA and A = 10~ 5 GeV 
for RP. The values of the coupling /3 are shown in the frames. 
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In Figure || we also plot w e ff, defined in accordance with eq. Q2.13 ). As expected from 
the required expansion rate, at large z it is w e ff = 1/3 . When w falls down and intersects 
w e ff, it also starts to decrease. In the inner frame of the Figure, the low— z behavior is 
magnified. Up to z ~ 1 the scale dependence of w and w e ff will be distinguishable only 
through very refined experiments. Above z ~ 1, however, the difference becomes more 
relevant. 

The interpolation ( |6.2| ) with Z\ = 11 and e = 1.9, then yields the behavior of the 
different cosmic component shown in Figure ||. Let us specifically outline that the final DE 
density can be tuned rather easily to different values; more specifically, a greater (smaller) 
Z\ yields a greater (smaller) today's DE density. On the contrary, DM cou final density is 
substantially insensitive to the parameter choice. 

Altogether, reproducing the observational densities requires a suitable set of model 
parameters, but no fine tuning is required: their tuning must be so precise as the parameter 
precision required. In the case of Figure 6, the present DE density is 3 times DM unc 
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Figure 5: Solid lines: DE state parameter during the kinetic-potential transition. Dashed lines: 
Effective state parameter, if hypothetical data are considered by supposing DE not to be coupled. 
The large j3 value balances the small DM cou density, so that the w scale dependence is not so 
different, e.g., from a RP case with f} = 0.1 and all DM coupled (Figure 4). The inner panel show 
the low— z behavior, outlining that w and w e ff are not so different up to z ~ 1, again because of 
the small DM cou density parameter. 



(including also baryons), while the final contribution of DM Cl 
the per mil range. 
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7. Astrophysical context: flat galaxy cores and small halo deficit 

If we then compare the cosmological picture proposed here with more standard scenar- 
ios, where DE is however a scalar field, the main difference is that here the (ft field has 
substantially contributed to the cosmic budget since ever. 

On the contrary, multiple DM components have been considered by several authors 
(see, e.g., [ 17| , [ 18 ] ) , for various reasons. In this Section we shall debate the option 
that DM unc is WDM, instead of CDM. In association with the presence of a second DM 
component, and owing to the peculiar features of DM cou , this option seems particularly 
appealing. 

Let us then first remind that, if a DM component is coupled to DE, the effective gravity 
it feels is modified. For baryonic matter, we have sophisticated tests allowing us to state 



that, on terrestrial and planetary scales, no coupling with DE exists [19]. No such test 
can be extended to DM and this is why quite a few options for DM-DE coupling were 
considered in the literature, and tested against cosmic data. 
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Figure 6: Density evolutions for a model as in Figure 2, with the DE state equation in Figure 5. 



Observations, however, so nicely fitting most ACDM predictions, put stringent limits to 
the coupling f3 [20]. They can be eased when DM-DE coupling is considered in association 
with non-vanishing neutrino masses and Mildly Mixed Coupled cosmologies were 
found to fit data slightly better than ACDM. The likelihood improvement, however, is 
not statistically significant and, even within this approach, a coupling range including 
j3 > 0.886 is excluded. If two DM components exist, however, we are in a fully different 
context, that we may tentatively exploit to seek a solution to some inconsistencies of the 
ACDM model, on sub-galactic scales, put in evidence by N-body simulations, namely if 
DM is assumed to be "cold" . 



A first difficulty concerns the amount of substructure in Milky Way sized haloes |23| . 
Models involving CDM overpredict their abundance by approximately one order of magni- 
tude. A second issue concerns the density profiles of CDM haloes in simulations, exhibiting 
a cuspy behavior 24, 25 1, while the density profiles inferred from rotation curves suggest a 
core like structure |2q|. A third issue concerns dwarf galaxies in large voids: recent studies 



1 27] re-emphasized that they are overabundant. 

It is known that replacing CDM with a "warmer" DM component, as a thermal relic of 
particles whose mass is ~ 2-3 keV, yields better predictions. What is essential, however, is 
the streaming length of such component. Accordingly, it can also be replaced by particles 
of different mass, with a non-thermal distribution, but similar average velocity. However, 
there is a number of "thermal" candidates for such warm dark matter (WDM); among 
them, a sterile neutrino and a gravitino [pi] find a reasonable motivation in particle theory 



The long streaming length of such particles causes a strong suppression of the power 
spectrum on galactic and sub-galactic scales |3(]] and solves several above problems. In 
particular, the profiles of WDM haloes, similar to CDM haloes in the outer regions, flatten 
towards a constant value in the inner regions, as predicted in |H]] and found in simulations 
& 

However, the core size found is 30-50 pc, while the observed cores in dwarf galaxies are 
around the 1000 pc scale [33]. A dwarf galaxy core in this scale range would be produced by 
higher velocity particles, as those belonging to a thermal distribution if their mass is < 0.1- 
0.3 keV. Increasing the velocity, however, yields a greater streaming length, exceeding the 
size of these very dwarf galaxies, in the first place [63]. 

In view of these difficulties, the idea that WDM is accompanied by a smaller amount 
of CDM has already been put forward [16]. The WDM particle velocities could then 
be greater, while a low-mass population is however produced by CDM clustering. This 
suggestion was been put forward quite independently of any particle or cosmic model. In 
particular, assuming ad hoc a twofold dark matter component does not ease coincidence 
problems. 

It is then clear that the model discussed in this paper could been adapted to meet 
the above requirement. DM unc would be a kind of WDM, made of high-speed particles. 
DM cou , then, would be responsible to create condensation sites on scales smaller than the 
DM unc streaming length, after its derelativisation. Its role is similar to the CDM role in 
ACDM models, after recombination, when CDM fluctuations cause baryon accretion on 
scales where primary baryon fluctuations had been erased during recombination. 

At large scale we then expect a standard primeval fluctuation spectrum, suitably bal- 
anced between DM cou and DM unc , although today's DM cou contribution could be not 
so significant. Below the streaming length of DM unc , however, only DM cou fluctuation 
initially remain, to create the seeds for the low scale fluctuation spectrum, after DM unc 
derelativization. The amplitude of the overall DM spectrum can then be expected to be a 
few times smaller below the WDM (DM unc ) streaming length. 

These qualitative considerations require a detailed quantitative confirm. It is not un- 
reasonable, however, that "secondary" DM unc fluctuations, although still yielding some 
low-mass structure, generate less low-mass haloes than a standard CDM spectrum. The 
evolution of the fluctuations in DM cou needs however a direct inspection. It is known 
that, at the Newtonian level, its gravitational self-interaction force is enhanced by a factor 
1 +4/3 2 /3 (~ 9 for j3 ~ 2.5). Therefore, after entering the non-linear regime, DM cou fluc- 
tuations, more rapidly than fluctuations in other components, could evolve into collapsed 
objects; their expected features are not easily predictable without a detailed analysis and 
should be then compared with observations. 

Analogous comments can be made for the size of the cores in low mass galaxies. Be- 
ing mostly made of low-mass WDM they can be as large as required, while the galaxy 
population does exist thanks to the DM cou spectral seeds. 

To put these expectations in a quantitative form we need to study fluctuation evolution 
in detail. 



8. Discussion 



In this work we considered coupled DE theories, when the coupling constant (3 is large. 
The first finding is then that a dual component, made of non-relativistic particles and a 
scalar field, can be in equilibrium with the radiative components in the radiative era. The 
density parameters of the dual components have then a fixed ratio 

n c /n d = 2, while (n c + n d )/3 2 = 3/4 (8.1) 

(if the scalar field energy is prevalently kinetic) and such density parameters keep constant, 
as both dual components dilute oc a -4 as the Universe expands. We dubbed them DM cou 
and DE although, at this stage, there is no evidence of the latter being related to observa- 
tional DE. Another important finding is that the dual component is stable: if we set initial 
conditions violating (|3.1|), the densities of DM COM and DE change and the condition ( |S.1| ) 
is restored. 

The dual component gives place to a natural form of DR. Here we find that 

AN °» = 3^075 ' (8 - 2 » 

so that j3 = 2.5 yields AN Q ff ~ 1. All plots in this paper are given for this case. 

When the expansion regime ceases to be radiative, as the density of a different non- 
relativistic component (oc a~ 3 ) overcomes the radiative component, also the densities of 
DM con and DE start to increase. This further DM component is dubbed DM UT1C . DM co1i 
and DE densities, although increasing, however keep well below DM unc , unless the energy 
of the scalar field shifts from kinetic to potential. 

Let us then recall again that, when we try to fit background data to a model where DE 
is coupled to the whole DM and self-interacts through a standard potential (e.g., Ratra- 
Peebles or SUGRA), we find that the 4> field energy shifts from kinetic to potential at 
a redshift Z\ ~ 10-20. This is the range where the transition must occur, also in this 
case. The example shown in the Figures is for z\ = 11, but similar proportions of DE are 
obtained also by slightly shifting z\ and the parameter e, simultaneously. Large shifts of 
such parameters are however not allowed and the epoch of the field transition from kinetic 
to potential is well constrained, quite independently of the freedom we still have to modify 
its detailed dynamics. 

Let us finally stress that we describe a stationary high-z situation, holding since the 
decoupling of the DM COM component from the other particles and, possibly, even before this 
stage. The dynamics discussed here could originate from a phenomenological Lagrangian 
coupling/mass term 

C = Tdnbip or C = , (8.3) 

m p 

holding since then. Here t/j describes spinor particles, which are DM cou . Should the validity 
of this Lagrangian extend until the end of the inflationary era, one wonders whether any 
relation exist between the DE (/>-field and the scalar field responsible for the inflationary 
process itself. 
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